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09/609,219, filed on June 30, 2000, now pending. 

FIELD AND BACKGROUND OF THE INVENTION: 

The present invention relates to the field of genomics, and more 
particularly, to a method and system for evaluating the electrical 
conductivity of regions of DNA. The well known Watson & Crick model of 
the double helix describes the DNA molecule essentially as an elongated 
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Stack of aromatic nitrogenous base pairs wrapped with a ribbon of a 
negatively charged sugar phosphate polymer pair. Carbon atoms in the ring 
structures of the base pairs occur in SP2 hybridization and have n orbitals 
perpendicular to the ring planes. The intimate packing of adjacent base pairs 
5 in the core stack suggests a high degree of overlap between their % electron 
clouds. Since shortly after the Watson & Crick model was first proposed 
(Watson JD and Crick FHC (1953) Nature 171:737), scientists have begun 
entertaining the idea that the structure of the DNA molecule suggests that it 

ti 

may possess electric conductivity (Bumel et al. (1969) Ann N Y Acad Sci 
10 158: 191-209; Holmlin RE (1997) Angew Chem Int Ed Engl 36, 
2714-2730). This idea has been extensively developed and tested in recent 
years, both experimentally (Aich P et al. (1999) J Mol Biol 294: 477-85; 
al-Kazwini AT et. al. (1994) Radiat Res 138: 307-1 1 ; Arkin MR et al. (1996) 
Science 273: 475-80; Arkin MR et al. (1997) Chem Biol 4: 389-400; Beach 
15 C et al. (1994) Radiat Res 137: 385-93; Dandliker PJ et al. (1997) Science 
275: 1465-8; Dandliker PJ et al. (1998) Biochemistry 37: 6491-502; Das A 
et al. (1999) Chem Biol 6: 461-71; Debije MG and Bemhard WA (1999) 
Radiat Res 152: 583-589; Debije MG et al. (1999) Angew Chem Int Ed Engl 
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38: 2752-2756; Douki T and Cadet J (1999) Int J Radiat Biol 75: 571-81; 
Femandez-Saiz M et al. (1999) Photochem Photobiol 70: 847-52; Fink HW 
and Schonenberger C (1999) Nature 398: 407-10; Fuciarelli AF et al. (1994) 
Int J Radiat Biol 65: 409-18; Gasper SM (1997) J Am Chem Soc 119: 
5 12762-12771; Hall DB (1997) J Am Chem Soc 119: 5045-5046; Hall DB et 
al. (1996) Nature 382: 731-5; Hall DB et al. (1998) Biochemistry 37: 
15933-40; Henderson PT et al. (1999) Proc Natl Acad Sci U S A 96: 8353-8; 
Holmlin RE (1997) Angew Chem Int Ed Engl 36: 2714-2730; Hyun KM et 
al. (1997) Biochim Biophys Acta 1334: 312-6; Jovanovic SV and Simic MG 
10 (1989) Biochim Biophys Acta 1008: 39-44; Kelley SO and Barton JK 
(1998) Chem Biol 5: 413-25; Kelley SO and Barton JK (1999) Science 283: 
375-81; Kelley SO et al. (1999 Nucleic Acids Res 27: 4830-4837; Kielkopf 
CL et al. (2000) Nat Struct Biol 7: 117-21; Lewis FD et al. (1997) Science 
277: 673-6; Meggers E et al. (2000) Chemistry 6: 485-92; Meggers E. (1998) 
15 J Am Chem Soc 120: 12951-12955; Murphy CJ et al. (1994) Proc Natl Acad 
Sci U S A 91: 5315-9; Murphy CJ et al. (1993) Science 262: 1025-9; Napier 
ME et al. (1997) Bioconjug Chem 8: 906-13; Ninomiya K et al. (1999) 
Nucleic Acids Symp Ser 255-6; Nunez ME and Barton JK (2000) Curr Opin 
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Chem Biol 4: 199-206; Nunez ME et al. (1999) Chem Biol 6: 85-97; Nunez 
ME et al. (2000) Biochemistry 39: 6190-6199; Ozaki H and McLaughlin LW 
(1992) Nucleic Acids Symp Ser 67-8; Porath D et al. (2000) Nature 403: 
635-8; Ratner M (1999) Nature 397: 480-1; Razskazovskiy Y et al. (2000) 
5 Radiat Res 153: 436-41; Ropp PA and Thorp HH (1999) Chem Biol 6: 
599-605; Schuster GB (2000) Acc Chem Res 33: 253-260; Thompson M and 
Woodbury NW (2000) Biochemistry 39: 4327-4338; Wagenknecht HA et al. 
(2000) Biochemistry 39: 5483-5491; Wan C et al. (1999) Proc Natl Acad Sci 
U S A 96: 6014-9; Wolf P et al. (1993) Int J Radiat Biol 64: 7-18; Xu DG 
10 and Nordlund TM (2000) Biophys J 78: 1042-58) and theoretically (Bixon 
M et al. (1999) Proc Natl Acad Sci USA 96: 11713-6; Buhks E and Jortner J 
(1980) FEBS Lett 109:117-20; Chojnacki H and Laskowski Z (1985) J 
Biomol Struct Dyn 2:759-65; Conwell EM and Rakhmanova SV (2000) Proc 
Natl Acad Sci USA 97:4556-60; Fonseca Guerra C and Bickelhaupt FM 
15 (1999) Angew Chem Int Ed Engl 38:2942-5; Grinstaff MW (1999) Angew 
Chem Int Ed Engl 38:3629-35; Jortner J et al. (1998) Proc Natl Acad Sci 
USA 95:12759-65; Meade TJ (1996) Met Ions Biol Syst 32:453-78; Nunez 
ME and Barton JK (2000) Curr Opin Chem Biol 4:199-206 ; Onfelt B et al. 
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(2000) Proc Natl Acad Sci USA 97: 5708-13; Schuster GB (2000) Acc 
Chem Res 33:253-60; Steenken S (1997) Biol Chem 378:1293-7; Steenken S 
(1992) Free Radic Res Commun 16:349-79; Takeda K (1995) Math Biosci 
130:183-202; Wan C et al. (1999) Proc Natl Acad Sci USA 96:6014-9.) 

5 Recent work carried out in a number of independent laboratories, 

using various biophysical methodological approaches, demonstrated that 
long-range migration of charge through the double helix occurs (Aich P et al. 
(1999) J Mol Biol 294: 477-85; al-Kazwini AT et. al. (1994) Radiat Res 138: 
307-1 1 ; Arkin MR et al. (1996) Science 273: 475-80; Arkin MR et al. (1997) 

10 Chem Biol 4: 389-400; Beach C et al. (1994) Radiat Res 137: 385-93; Clery 
D (1995) Science 267:1270; Dandliker PJ et al. (1997) Science 275: 1465-8; 
Dandliker PJ et al. (1998) Biochemistry 37: 6491-502; Das A et al. (1999) 
Chem Biol 6: 461-71; Dunn DA et. al., (1992) Biochemistry 31: 11620-5; 
Fink HW and Schonenberger C (1999) Nature 398: 407-10; Hall DB et al. 

15 (1996) Nature 382: 731-5; Hall DB et al. (1998) Biochemistry 37: 15933-40; 
Henderson PT et al. (1999) Proc Natl Acad Sci U S A 96: 8353-8; Holmlin 
RE (1997) Angew Chem Int Ed Engl 36: 2714-2730; Jovanovic SV and 
Simic MG (1989) Biochim Biophys Acta 1008: 39-44; Murphy CJ et al. 
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(1994) Proc Natl Acad Sci U S A 91: 5315-9; Murphy CJ et al. (1993) 
Science 262: 1025-9; Nunez ME and Barton JK (2000) Curr Opin Chem 
Biol 4: 199-206; Nunez ME et al. (2000) Biochemistry 39: 6190-6199; 
Porath D et al. (2000) Nature 403: 635-8; Wagenknecht HA et al. (2000) 
5 Biochemistry 39: 5483-5491; Wan C et al. (1999) Proc Natl Acad Sci U S A 
96: 6014-9) and that it is sequence and n stacking sensitive (Arkin MR et al. 

(1996) Science 273: 475-80; Fuciarelli AF et al. (1994) Int J Radiat Biol 65: 
409-18; Hall DB (1997) J Am Chem Soc 119: 5045-5046; Kelley SO and 
Barton JK (1998) Chem Biol 5: 413-25; Kelley SO and Barton JK (1999) 

10 Science 283: 375-81; ; Meggers E et al. (2000) Chemistry 6: 485-92; 
Meggers E. (1998) J Am Chem Soc 120: 12951-12955; Napier ME et al. 

(1997) Bioconjug Chem 8: 906-13; Nunez ME and Barton JK (2000) Curr 
Opin Chem Biol 4: 199-206; Nunez ME et al. (1999) Chem Biol 6: 85-97; 
Ropp PA and Thorp HH (1999) Chem Biol 6: 599-605; Schuster GB (2000) 

15 Acc Chem Res 33: 253-260; XU DG and Nordlund TM (2000) Biophys J 
78: 1042-58.) While most of the direct evidence supporting DNA charge 
migration was obtained in highly artificial experimental systems, the impact 
of these findings on biological research may be immense because they raise 



the possibility that the double helix may have a capacity to conduct charge 
also under normal physiological conditions (Blank M and Goodman R 
(1997) Bioelectromagnetics 18:111-5; Blank M and Goodman R (1999) J 
Cell Biochem 75: 369-74; Buhks E and Jortner J (1980) FEBS Lett 
109:117-20; Cullis PM et. al. (1987) Nature 330:773-4.) 

DNA is the universal carrier of genetic information. If it is true that 
charge transfer through the double helix exists under normal physiological 
conditions in cells of living organisms, it could have a profound impact on 
the understanding of the function of genetic material and processes which 
take place directly on it (such as transcription and replication). It would not 
be impossible that charge transfer in the double helix is utilized by the cell as 
a component of the control mechanisms that govern gene expression. 

In light of the results demonstrating DNA charge migration in vitro, 
prominent researchers in the field have expressed the opinion that attempts 
should be directed at studying the biological implications of the phenomenon 
(Blank M and Goodman R (1997) Bioelectromagnetics 18:111-5; Blank M 
and Goodman R (1999) J Cell Biochem 75: 369-74; Grinstaff MW (1999) 
Angew Chem Int Ed Engl 38:3629-35; Holmlin RE (1997) Angew Chem Int 
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Ed Engl 36: 2714-2730; Lin H et. al. (1998) J Cell Biochem 70: 297-303; 
Nunez ME and Barton JK (2000) Curr Opin Chem Biol 4: 199-206; Nunez 
ME et al. (1999) Chem Biol 6: 85-97.) If charge migration serves a 
biological function, then traces of the phenomenon may be predicted to have 

5 been registered into the genome by evolution. To date, there has been no 
publication of any experimental approach for rigorous testing of hypotheses 
on the putative biological effects of DNA charge migration, nor has there 
been any publication of a method for determining the electrical conductivity 
properties of a selected DNA segment. 

10 Biophysical evidence (see references hereinabove) suggests that 

charge migration is sequence sensitive. If this is true then the charge 
conductivity qualities of a given region in a chromosome are dictated by the 
sequence of bases. If a particular region in genomic DNA functions as a 
modulator of charge conductivity, its structure can be expected to reflect that 

15 fact. Non-coding regions of genomic DNA do not code for proteins but they 
still code for their own physicochemical characteristics. In DNA 
'structure = sequence' and structural characteristics of the molecule can 
therefore be deduced from sequence analysis. Thus, depending on its specific 
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sequence, the electrical conductivity properties of a region of DNA might be 
predicted. There is no prior art suggesting that, or how, this might be 
accomplished. While prior art methods and systems for analysis of DNA 
sequence exist, there are no prior art systems that provide a method or means 

5 for analysis of a DNA sequence for determination of its electrical 
conductivity properties. 

For example, prior art systems (Brendel V and Trifonov EN (1984) 
Nucleic Acids Res 12:4411-27; Cox R and Mirkin SM (1997) Proc Natl 
Acad Sci USA 94:5237-42; Day GR and Blake RD (1982) Nucleic Acids 

10 Res 10:8323-39; Karlin S (1986) Proc Natl Acad Sci USA 83:6915-9; Karlin 
S and Brendel V (1992) Science 257:39-49; Karlin S et. al. (1992) Nucleic 
Acids Res 20:1363-70; Karlin S and Cardon LR (1994) Annu Rev Microbiol 
48:619-54; Karlin S and Ghandour G (1985) J Mol Evol 22:195-208; Karlin 
S et. al. (1983) Proc Natl Acad Sci USA 80:5660-4; Karlin S et. al. (1988) 

15 Comput Appl Biosci 4:41-51; Karlin S et. al. (1988) Proc Natl Acad Sci 
USA 85:841-5; Kom LJ et. al. (1977) Proc Natl Acad Sci USA 74:4401-5; 
Mount DW and Conrad B (1984) Nucleic Acids Res 12:811-7; Munroe SH 
(1983) Nucleic Acids Res 11:8891-900; Nussinov R (1991) DNA Seq 
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2:69-79; Pivek et.al. (1985) Folia Biol (Praha) 31:213-34; Schroth GP and 
Ho PS (1995) Nucleic Acids Res 23:1977-83) have been described which 
utilize algorithms to analyze DNA sequences for restriction sites, 
homologies to other sequences, or homologous elements such as repetitive 
elements and complete dyad symmetries and palindromes, in order to 
compare genomes of different species, and to identify potential cruciform 
and hairpin loop secondary structures. These various computer algorithms, 
based on different methods of operation, have been devised to locate dyad 
symmetrical elements in sequence and determine their frequency and 
distribution in order to predict secondary structure of the nucleic acid 
molecules. However, all these systems suffer from limitations, such as that 
they are designed to identify symmetry and not asymmetry (generally 
functioning by attempting to locate two complementary and reversed 
occurrences of the same nucleotide sequence in close proximity on the same 
strand.) They were not designed to or able to quantitate the degree of 
asymmetry, even though some are tolerant to minor degrees of mismatch. 
None provides a means for analysis of a DNA sequence for determination of 
its electrical conductivity properties. 
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Most of the work on fimctional characterization of genomic DNA 
sequence (locating functional elements such as promoters, enhancers, splice 
sites, and origins of replication, for example) is empirical. Most of the data 
in the field is obtained through elaborate and time consuming in- vivo and in 
vitro assays. A need exists for computer based tools, which enable scientists 
to accurately predict functional elements in DNA sequence by means of 
sequence analysis. Tools of this kind can reduce the amount of costly and 
time consuming 'wet' experiments, and enhance the efficiency of the process 
of understanding the complexity of the genome. Effective tools for 
assembling short fragments of expressed sequence into putative exons are in 
widespread commercial and academic use, providing an example of the need 
and demand for such tools in general. A particular demand exists for new 
and improved tools that would enable prediction of the location and strength 
of elements involved in gene expression, such as promoters and enhancers 
with enhanced accuracy and precision. Such improved tools could reduce the 
amount of labor required by current methodologies for determination of the 
location of transcription start sites and for evaluation of the strength of 
promoters and enhancers. It could also be used for the study and diagnosis of 
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various genetic diseases and conditions, where these are caused by mutations 
in the regulatory elements of genes rather than by mutations in protein 
coding sequences. 

Thus there is an unmet need for, and it would be highly advantageous 
to have, a method for identifying DNA sequences that are good candidates to 
effectively propagate charge transfer and serve as an electrical conductor. 
There is further an unmet need for, and it would be highly advantageous to 
have, computer based tools to enable scientists to accurately predict 
fiinctional elements, such as transcription related functional elements, in 
DNA sequence by means of sequence analysis. 



SUMMARY OF THE INVENTION: 

According to the present invention there is provided methods and a 
system to determine the electrical conductivity properties of a DNA 
sequence. Further, there is provided a method and a system for identifying 
functional elements in a DNA sequence. 
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According to one aspect of the present invention there is provided a 
method for determining a measure of electrical conductivity of a defined 
DNA sequence, which comprises the step of calculating the degree of 
asymmetry of the defined DNA sequence. 

According to another aspect of the present invention there is provided 
a system for determining a measure of electrical conductivity of a defined 
DNA sequence, which comprises a computer containing within a memory 
device thereof, an algorithm which is capable of calculating the degree of 
asymmetry of the defined DNA sequence. 

According to yet another aspect of the present invention there is 
provided a method for the evaluation of electrical conductivity of a defined 
DNA sequence which comprises the step of providing instructions on a 
computer readable medium for a calculation of the degree of asymmetry of 
the defined DNA sequence. 

According to still another aspect of the present invention there is 
provided a method for identifying functional elements in a DNA sequence, 
the method including the steps of: (a) calculating at least one set of dyad pair 
type frequencies within a portion of the DNA sequence equal in size to two 



times a window size, around at least one potential axis of dyad symmetry in 
the DNA sequence, and (b) based on the at least one set of dyad pair type 
frequencies, identifying regions of the DNA sequence containing the 
functional elements. 

5 According to an additional aspect of the present invention there is 

provided a method of identifying transcription related functional elements, 
including the steps of: (a) calculating at least one set of dyad pair type 
frequencies within a portion of the DNA sequence equal in size to two times 
a window size, around at least one potential axis of dyad symmetry in the 

10 DNA sequence, and (b) based on the at least one set of dyad pair type 
frequencies, identifying regions of the DNA sequence containing the 
functional elements. 

According to yet an additional aspect of the present invention there is 
provided a system for identifying functional elements in a DNA sequence, 

15 the system including: (a) a software module including a plurality of 
instructions for calculating at least one set of dyad pair type frequencies 
within a portion of the DNA sequence equal in size to two times a window 
size, around at least one potential axis of dyad symmetry in the DNA 
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sequence, (b) a memory for storing the instructions, and, (c) a processor for 
executing the instructions. 

According to still an additional aspect of the present invention there is 
provided a computer readable storage medium having computer readable 

5 code embodied on the computer readable storage medium, the computer 
readable code for identifying functional elements in a DNA sequence, the 
computer readable code comprising: program code including a plurality of 
instructions for calculating at least one set of dyad pair type frequencies 
within a portion of the DNA sequence equal in size to two times a window 

10 size, around at least one potential axis of dyad symmetry in the DNA 
sequence 

According to still an additional aspect of the present invention there is 
provided a method for determining electrical conductivity properties of a 
DNA sequence, the method comprising the steps of: (a) calculating at least 
15 one set of dyad pair type frequencies within a portion of the DNA sequence 
equal in size to two times a window size, around at least one potential axis of 
dyad symmetry in the DNA sequence, and, (b) based on the at least one set 
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of dyad pair type frequencies, determining the electrical conductivity 
properties of the DNA sequence. 

According to still an additional aspect of the present invention there is 
provided a system for determining electrical conductivity properties of a 
DNA sequence, the system including: (a) a software module including a 
plurality of instructions for calculating at least one set of dyad pair type 
frequencies within a portion of the DNA sequence equal in size to two times 
a window size, around at least one potential axis of dyad symmetry in the 
DNA sequence, (b) a memory for storing the instructions, and, (c) a 
processor for executing the instructions. 

According to yet an additional aspect of the present invention there is 
provided a computer readable storage medium having computer readable 
code embodied on the computer readable storage medium, the computer 
readable code for determining electrical conductivity properties of a DNA 
sequence, the computer readable code comprising: program code including a 
plurality of instructions for calculating at least one set of dyad pair type 
frequencies within a portion of the DNA sequence equal in size to two times 
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a window size, around at least one potential axis of dyad symmetry in the 
DNA sequence. 

According to further features in preferred embodiments of the 
invention described hereinbelow, the calculation of degree of asymmetry is 
accomplished by calculating at least one polarity value around at least one 
potential axis of dyad symmetry in the DNA sequence, where the polarity 
value is a unit-less number defined as (1-[S/W]), where S represents a 
number of dyad- symmetrical bases and W represents a window size. 

According to still further features in preferred embodiments of the 
present invention, the at least one polarity value is an ordered series of 
polarity values iteratively calculated for each potential axis in the DNA 
sequence. 

According to still further features in preferred embodiments of the 
present invention, the series of polarity values is plotted graphically, 
whereby extended regions in the DNA sequence that possess values of 
polarity which deviate from expected polarity values of a random sequence 
may be identified. 
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According to still further features in preferred embodiments of the 
present invention, the series of polarity values is subjected to statistical 
analysis, whereby extended regions in the DNA sequence that possess values 
of polarity which deviate from expected polarity values of a random 
sequence may be identified. 

According to still further features in the described preferred 
embodiments of the present invention, the DNA sequence length is in the 
range of 2 bases to 3 X 10^ bases. 

According to still further features in preferred embodiments of the 
invention described below, the window size is an independent variable, with 
values ranging from 1 to a value equal to that of the largest whole integer 
smaller than one half the length of the DNA sequence, to be designated prior 
to calculation. 

According to still further features in preferred embodiments of the 
invention described below, at least one of the at least one set of dyad pair 
type frequencies is an ordered array of the dyad pair frequencies. 

According to still further features in preferred embodiments of the 
invention described below, the calculating of at least one set of dyad pair 
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type frequencies is effected iteratively for each of the at least one potential 
axis of dyad symmetry in the DNA sequence. 

According to still further features in preferred embodiments of the 
invention described below, the step of identifying regions of the DNA 
sequence containing the functional elements is effected by steps including 
subjecting the at least one set of dyad pair type frequencies to statistical 
analysis, whereby at least one region in the DNA sequence is identified that 
possesses at least one statistical value that indicates that an observed at least 
one set of dyad pair type frequencies deviates from an expected at least one 
set of dyad pair type frequencies. 

According to still fiirther features in preferred embodiments of the 
invention described below, the at least one statistical value is chosen from 
the group consisting of residuals of the dyad pair type frequencies, 
chi-square values, and likelihood ratios. 

According to still further features in preferred embodiments of the 
invention described below, the statistical analysis includes plotting said 
statistical values. 
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According to still further features in preferred embodiments of the 
invention described below, the window size is an independent variable, 
having a value of at least 1 and at most one half a length of the DNA 
sequence. 

According to still further features in preferred embodiments of the 
invention described below, the calculating of the at least one set of dyad pair 
type frequencies is performed on at least one dyad pair wherein both 
nucleotides of the at least one dyad pair are located on a single strand of the 
DNA sequence. 

According to still further features in preferred embodiments of the 
invention described below, the calculating of the at least one set of dyad pair 
frequencies is performed on at least one dyad pair wherein both nucleotides 
of the at least one dyad pair are located on complementary strands of the 
DNA sequence. 

The method and system of the present invention can be used as a 
research instrument to enable the search for evidence supporting the theory 
of DNA charge migration in DNA sequence. If charge migration serves a 
biological function, then traces of the phenomenon may be registered into 
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the base sequence by evolution. As will be shown below, examination of 
specific DNA sequences, such as the genome of the human 
immunodeficiency virus, HIV2, with this method and system has been 
accomplished yielding the unexpected result that several apparent extended 

5 regions of increased polarity are easily identified. 

When applied to genomic DNA sequence of various organisms, 
particularly human, the method can locate distinct, recognizable elements in 
DNA. When the location of these distinct elements is superimposed on the 
functional map of genes, it is evident that a significant degree of overlap 

10 exists. One striking example is the element that is found to reside on the 
transcription initiation point of many human genes. The system and method 
of the present invention, by detecting this element, can accurately predict the 
location of promoters of genes. It appears that not only the location but also 
the strength of a promoter can be predicted. Some genetic diseases and 

15 conditions, such as Fragile-X syndrome, for example, are caused by 
mutations of the control regions of genes rather than of their protein coding 
regions. When affected and non-affected alleles of the FMRl gene (known 
to be responsible for Fragile-X syndrome) are analysed using the present 
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invention, a clearly visible difference is found between the transcription 
initiation site elements of the affected and non-affected alleles. As many 
genetic conditions may be caused by mutations in control regions, the system 
and method of the present invention will therefore contribute to a better 
understanding of the control of gene expression in both normal and 
pathological situations. 

The present invention thus successfully addresses the shortcomings of 
the presently known configurations by providing a method and a system for 
the analysis of a defined nucleotide sequence to calculate the degree of 
asymmetry in that sequence in order to determine the electrical conductivity 
properties of a DNA sequence. The present invention further provides a 
method and a system for identifying functional elements within a DNA 
sequence. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is herein described, by way of example only, 
with reference to the accompanying drawings. With specific reference now 
to the drawings in detail, it is stressed that the particulars shown are by way 
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of example and for purposes of illustrative discussion of the preferred 
embodiments of the present invention only, and are presented in the cause of 
providing what is believed to be the most useful and readily understood 
description of the principles and conceptual aspects of the invention. In this 
regard, no attempt is made to show structural details of the invention in more 
detail than is necessary for a fundamental understanding of the invention, the 
description taken with the drawings making apparent to those skilled in the 
art how the several forms of the invention may be embodied in practice. 
In the drawings: 

FIG. lA is a drawing that schematically illustrates a fragment 
of duplex DNA, 16 base pairs long (an example of an input sequence); along 
the top of the figure, potential axes of symmetry are indicated; 

FIG. IB is a flow diagram indicating the major steps in the 
analysis of the input sequence; for illustrative purposes only, the window 
size is set here to five nucleotides; dyad- symmetrical bases are indicated by 
bold typeface; each iterative step is representative of the application of the 
formula to calculate polarity at a different potential axis of symmetry; the 
input sequence and potential axes of symmetry are as indicated in figure 1 A; 
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FIG. IC is a table listing the polarity values for the 13 potential 
axes calculated in an example analysis of the input sequence in figure lA, 
using a window size of five, according to steps as illustrated in figure IB; 

FIG. ID is a graphic presentation of the output list of polarity 
values, taken from the example of figure IC; the expected value of 0.75 is 
also indicated; 

FIG. 2A is a graphic presentation of the output list of polarities 
from an analysis of the complete genome of HIV2 using an embodiment of 
the present invention; extended regions of increased polarity are indicated, 
these being regions of 500-600 nucleotides where values of polarity are 
concentrated which deviate from the 0.75 expected, and thus determined to 
be regions that would function effectively to propagate electron transfer and 
serve as an electrical conductor; 

FIG. 2B is a graphic presentation of the output list of an 
analysis of a randomly generated mock DNA sequence illustrating no 
significant extended deviation from the expected 0.75 value. 

FIG. 3 is a flow chart illustrating a preferred embodiment of 
the present invention; 
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FIG. 4 is a detailed flow chart of a computer algorithm further 
illustrating an example of a possible configuration of the present invention 
with iterative calculation of polarity values for muhiple potential axes of 
symmetry for a DNA sequence; 

FIG. 5 A is a table illustrating the 16 dyad pair types; 

FIG. 5B is a drawing that schematically illustrates a fragment 
of duplex DNA, 8 base pairs long (an example of an input sequence); a 
potential axis of symmetry at the center of the fragment and a window size 
of 4 are indicated; 

Fig 6 is a flow chart illustrating an alternate preferred 
embodiment of the present invention; 

Fig. 7 is a detailed flow chart of a computer algorithm further 
illustrating an example of a possible configuration of a preferred 
embodiment of the present invention with iterative calculation of dyad pair 
type frequencies for multiple potential axes of symmetry for a DNA 
sequence; 
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Fig. 8 shows output plots of dyad pair type frequency analyses 
of nine genomic sequence fragments and a control fragment of computer 
generated random sequence; 

Fig. 9 shows a dyad pair type frequency analysis of two DNA 
sequences, the FMRl fragment in FIG. 9A and the "FMR1+(CGG)333" 
fragment in FIG. 9B; and, 

Fig. 10 is a high level block diagram of a system for predicting 
the electrical conductivity properties and for identifying fimctional elements 
in a defined DNA sequence according to the present invention. 

DETAILED DESCRIPTION OF THE INVENTION: 

The present invention is of a method and system consisting of a 
computer algorithm which can be used to determine a measure of electrical 
conductivity of a defined DNA sequence. Specifically, the present invention 
can be used to calculate the defined DNA sequence's degree of asymmetry 
over an extended length. Furthermore, the present invention is used to 
identify functional elements, in particular transcription-related functional 
elements, within the DNA sequence. 
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The principles and operation of a method and system to determine a 
measure of electrical conductivity of a defined DNA sequence according to 
the present invention may be better understood with reference to the 
drawings and accompanying descriptions. For purposes of this specification 
and accompanying claims, the phrase, "dyad symmetry" is defined as the 
following: A DNA nucleotide sequence is said to show complete dyad 
symmetry when the base sequence at a particular position relative to an axis 
perpendicular to the DNA sequence on one strand of double-stranded DNA 
is identical to the base sequence on the complementary strand at a position 
equidistant from the axis, although in opposite orientation (that is, reading 
left to right on the upper strand for example, and right to left on the 
complementary lower strand). The degree to which the base sequence at a 
particular position relative to an axis perpendicular to the major longitudinal 
axis of DNA molecule on one strand of double-stranded DNA is identical to 
the base sequence on the complementary strand at a position equidistant 
from the axis indicates the degree of symmetry of that sequence if they are 
less than completely identical. Two bases are said to be dyad symmetric 
when the two bases, at the same position (distance) relative to an axis 
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perpendicular to the major longitudinal axis of DNA molecule, but located 
on opposite strands of double stranded DNA, are identical. For purposes of 
this specification and accompanying claims, in certain configurations of the 
present invention two bases may also be considered dyad symmetric (that is 
there is dyad symmetry present) when the two bases, at the same position 
(distance) relative to an axis perpendicular to the major longitudinal axis of 
DNA molecule, but located on opposite strands of double stranded DNA, are 
not identical, but both belong to the same family of bases, that is, both are 
either purines or pyrimidines. 

For purposes of this specification and accompanying claims, the 
phrase, "axis of symmetry" is defined as an axis perpendicular to the major 
longitudinal axis of DNA molecule around which the nucleotide sequence 
can be analyzed to determine the degree to which the nucleotide sequence on 
one strand is identical to the base sequence on the complementary strand at a 
position equidistant from the axis, although in opposite orientation (that is, 
reading left to right on the upper strand for example, and right to left on the 
complementary lower strand). Because dyad symmetry may or may not be 
present around any given axis chosen, the axis may preferably be referred to 
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as a potential axis of dyad symmetry. For the purposes of this specification 
and the accompanying claims, the terms "axis of symmetry", "axis of dyad 
symmetry," "potential axis of symmetry," and "potential axis of dyad 
symmetry" shall be interpreted as meaning the same thing. For purposes of 
this specification and accompanying claims, the phrase, "window of 
symmetry" or "window size" is defined as the length in bases of the 
sequence being tested for identity at each side of any potential axis of 
symmetry. 

Before explaining at least one embodiment of the invention in detail, 
it is to be understood that the invention is not limited in its application to the 
details of construction and the arrangement of the components set forth in 
the following description or illustrated in the drawings (for example, any 
particular software programming language). The invention is capable of 
other embodiments or of being practiced or carried out in various ways. 
Also, it is to be understood that the phraseology and terminology employed 
herein is for the purpose of description and should not be regarded as 
limiting. 
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Any given fragment of double stranded DNA has two complementary 
5'^3' sequences, one for each strand. While in some cases (i.e., in perfect 
palindromes) these sequences may be identical, in the majority of 
circumstances they are different from each other (see figure lA). Comparing 
the 5'->3' sequence of a DNA fragment to the 5'— >3' sequence of the 
complementary strand of the same fragment is equivalent to comparing the 
two paths that a hypothetical test charge migrating through the fragment of 
DNA in either direction could take. Analogous to a diode, if a region of 
DNA has evolved to function as a charge conductivity modulation element, 
it is unlikely to exert its action in both directions equally and its sequence is 
therefore predicted to show a distinct directionality. Such directionality in a 
sequence can be revealed by a systematic analysis of polarity, that is, the 
extent of sequence asymmetry of the complementary strands over an 
extended length of base pairs. Regions of DNA with enhanced charge 
conductivity will be identified as extended regions with increased polarity as 
compared with expected. Extended regions with decreased polarity as 
compared to expected are also identified and are also predicted to possess 
unique charge conductivity properties, namely high resistance. 
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The input to the algorithm is a string of characters representing the 
order of nucleotide bases from a single strand of a molecule of DNA. 
Depending on the specific embodiment of the present invention, the output is 
a number or a series of numbers each representing the polarity value of one 
potential axis of dyad symmetry in the input sequence. Polarity is a unit-less 
number defined as (1-(S/W)), where S= number of dyad-symmetrical bases 
and W= window (see figure IB). A perfect palindrome has zero polarity at 
its central axis of dyad symmetry, and a homogenous stretch of DNA 
consisting on one strand of only one of the four bases (i.e. 
AAAAAAAAA. . .A) has a polarity value of one. 

The algorithm of a preferred embodiment calculates polarity by 
comparing the nucleotide sequence of a specified window size (number of 
base pairs) upstream of the tested potential axis of dyad symmetry, against 
the nucleotide sequence of an equal size downstream from this axis on the 
complementary DNA strand (see figure IB). In an alternative preferred 
embodiment described in greater detail hereinbelow, the nucleotide sequence 
of a specified window size upstream of the tested potential axis of dyad 
symmetry is compared against the nucleotide sequence of an equal size 
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downstream from this axis on the same DNA strand. A further feature is that 
the algorithm may perform this routine for each potential axis of dyad 
symmetry along the input sequence (see figures lA and B) and return an 
ordered list of polarity values for all of the axes tested (see figure IC). A 
margin, equal in size to the window of symmetry, must be excluded from 
analysis at each end of the input sequence. This is due to the fact that if an 
axis is chosen within this margin, the size of the window will exceed the 
number of bases present on one DNA strand, between the axis and the end of 
the input sequence. There exists one potential axis of dyad symmetry on 
each nucleotide base and one between every two consecutive bases (see 
figure lA). 

According to further features of the preferred embodiment, once the 
list of polarity values for all individual potential axes of dyad symmetry in 
the tested sequence is obtained, its content may be displayed in a graph (see 
figure IC, and figure 2). The graph presents the polarity value at each 
potential axis of dyad symmetry (or a moving average of groups of axes) 
along the tested sequence as the y coordinate. The abscissa (x) values of the 
graph are the axis numbers and can be readily associated with nucleotide 
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positions on the input sequence. The expected polarity value for a random 
sequence is 0.75, based on both theoretical calculation and experimental data 
with randomly generated sequence (see figure 2B). According to still further 
features of the preferred embodiment of the present invention, statistical 
analysis can be performed on the list of polarity values. Statistical analysis 
can be performed to calculate a probability ratio indicating the deviation of 
the observed polarity values from the expected. Standard statistical methods 
which will be familiar to those ordinarily skilled in the art may be used (see 
Brezinski DP (1975) Nature 253:128-30.) The specific statistical method to 
be used may be tailored to different configurations of the present invention. 
For example, variations in base composition in different organisms and in 
different regions of the genome (must) warrant the use of different statistical 
evaluations. 

Detailed description of the steps involved in a preferred embodiment 
of the present invention will be easiest to understand with reference to the 
numbered steps on the algorithmic flowcharts in figures 3 and 4. Figure 3 is 
a flow chart illustrating a specific embodiment of the present invention, with 
an example of the steps an algorithm for determination of the degree of 
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asymmetry of a defined DNA sequence could take, while the flow chart in 
figure 4 illustrates a further, even more specific example, of a preferred 
embodiment of the present invention, in the form of an algorithm 
implemented in PERL programming language. The variable names and 
functions indicated in bold in figure 4 are used by way of example and no 
details in these examples should be taken as limiting the application of this 
invention. The first step (1) is for a nucleotide sequence of a single strand of 
DNA (input sequence, $input_seq) of a desired length to be input. The 
sequence may be of any length from two bases to 3 X 10^ bases, preferably 
from 5,000 to 50,000, and most preferably firom 10,000 to 20,000. The 
second step (2) is the input of length of the desired window size ($win_sym). 
Window size (W as described hereinabove) may be any number from one to 
a value equal to that of the largest whole integer smaller than one half the 
length of the DNA sequence, preferably from 20 to 300 and most preferably 
from 80 to 100. Beginning then at 3, the algorithm starts at the first potential 
axis of symmetry (axis position = 2*the window size) and calculates and 
outputs a polarity value for that axis (4). In the detailed example of figure 4, 
the input sequence is converted to the two complementary sequence indexed 
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arrays: @trgt_fwd and @trgt_revcomp in the steps indicated as 3 using 
string $win_seq in the process of these steps. The algorithm tests all the pairs 
of isometric bases within the window around that axis for identity. In figure 
4, each base within that window is indexed using variable $i and the number 
of identical bases (S, as described hereinabove) is counted in variable 
$match_count. Polarity is calculated using the formula as described 
hereinabove, polarity = (1-(S/W)). In the algorithm of the example in figure 
4, polarity is recorded as the variable $asym_count and output to an indexed 
array, @axis_list. In the next steps 5 and 5', the algorithm advances to the 
next potential axis; variables $basefeed and $basefeed_comp are used to 
advance the axis and sequence in the example of figure 4. In steps 6 and 6', 
polarity value around the new axis is calculated and output and steps 5/5' 
and 6/6' are repeated iteratively up to and including axis position = 2*length 
of the input sequence -(2*window size). In step 7 the ordered list of polarity 
values around each potential axis of symmetry is output. Graphical, step 8, 
and statistical (step 9) analysis can be performed, allowing for identification 
of extended regions of increased polarity, step 10. For example, it can easily 
be seen in figure 2A that such regions are easily identifiable. Extended 
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regions with decreased polarity as compared to expected are also identified 
and are also predicted to possess unique charge conductivity properties, 
namely high resistance. 

An alternative preferred embodiment performs a more detailed 
analysis of dyad symmetry. The two bases that are situated at the same 
position (distance) relative to an axis perpendicular to the major longitudinal 
axis of a DNA molecule, but located on opposite strands of double stranded 
DNA, are referred to as a dyad pair. Each dyad pair is one of 16 possible 
permutations of bases, as illustrated in Figure 5 A. Each of these 16 
permutations is referred to as a dyad pair type (DPT). The 16 DPTs can be 
grouped into four groups: self dyad, self mirror, purine-pyrimidine dyad, and 
purine-pyrimidine mirror, as illustrated in figures 5A and 5B. Fig. 5B 
illustrates an 8 base pair fragment of DNA, a potential axis of symmetry at 
the center of the fragment and a window size of 4. The dyad pairs, from the 
axis outward, are examples of self dyad (^ - g)? self mirror (^ - c), 
purine-pyrimidine dyad (^ - a), and purine-pyrimidine mirror C - t) DPTs, 
respectively. Thus the self mirror group, for example, consists of the dyad 
pairs: - c (as seen in the second dyad pair in Fig. 5B), -j, - a, and -q. 
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In an alternate preferred embodiment of the present invention, illustrated 
schematically in the flow chart of Fig. 6, at each potential axis of symmetry, 
rather than calculating symmetry and polarity, the algorithm calculates the 
frequencies of each of the 16 possible DPTs of the sequence within a 
fragment of sequence equal to twice the size of a defined window of 
symmetry, relative to the central axis of that fragment. The sum of the four 
DPT frequencies in the self dyad group is the same as the symmetry measure 
("s") in the preferred embodiment described hereinabove. The sum of the 
frequencies of the self mirror, purine-pyrimidine dyad, and 
purine-pyrimidine mirror groups together is equivalent to the polarity 
measure (p) in the preferred embodiment described hereinabove. Thus, this 
preferred embodiment gives finer resolution than the analysis of the 
preferred embodiment described hereinabove and illustrated in figs 3 and 4. 
After determining the set of DPT frequencies, and statistical measures (as 
described hereinbelow), at the first potential axis of symmetry in the input 
DNA sequence, the algorithm advances to the next potential axis of 
symmetry, reiterates the calculation of DPT frequencies and associated 
statistical measures, and moves on until the end of the input sequence is 
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reached. This is done in a manner analogous to that described hereinabove 
for the preferred embodiment illustrated in figures 3 and 4. An ordered array 
of DPT frequencies and statistical measures around each potential axis of 
symmetry is output. 

As described hereinabove for the preferred embodiment illustrated in 
figures 3 and 4, in the preferred embodiment illustrated in Figs. 6 and 7, 
statistical and graphical analysis can be performed, allowing for 
identification of regions predicted to possess unique charge conductivity 
properties and regions predicted to encode functional elements. To assess the 

statistical significance of observed deviations from expected DPT 

2 

frequencies a chi-square (% ) statistic may be calculated as a non-limiting 

2 

example. The x value for each axis and window is calculated according to 
the formula: 



where / is the DPT as indicated in fig. 5A, Obi is the observed DPT 

frequency for that DPT and Exi is the expected DPT frequency for that DPT. 

2 

In calculating the x value, expected DPT frequencies are preferably 
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calculated based on a model in which a probability of 1/(4^) is assigned to 
each DPT. This probability model is based on the assumption of unbiased 
nucleotide composition, as expected for a random sequence. Thus, Ex = 
(1/16) X W. In other configurations, the expected frequencies can be based 

5 on a model using actual base composition, as counted in each window, or as 
counted for the entire fragment (two windows on either side of the axis 
combined), or it can be based on actual base composition, as counted for a 
particular chromosome, part of a chromosome, or the whole genome of a 
particular organism. As further non-limiting examples, a residual of the DPT 

10 frequency, where Residual,- = {Obi - Ex^, or a likelihood ratio indicating the 
deviation of the observed DPT frequencies from the expected also can be 
calculated. The % values, residuals, likelihood ratios and DPT frequencies 
can be graphically plotted against their axis position in the input sequence. A 
fragment of computer generated random sequence subjected to the same 

15 analysis serves as a negative control and helps to verify that the observations 

2 

are not an artifact of the analysis and that the x value threshold used is 
appropriate. Examples of such graphical plotting are given in Figs. 8 and 9, 
which are discussed in greater detail hereinbelow. As discussed hereinabove, 
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in some configurations, and in some calculations, various dyad pair types 
may be taken together, such as the 4 major groups as a non- limiting 
example. In some configurations, some of the statistical analysis of DPT 
frequency deviation is performed at the time of each set of DPT frequency 
calculations at each axis rather than following the calculation of all DPT 
frequencies. 

Figure 6 is a flow chart illustrating a specific preferred embodiment of 
the present invention, with an example of the steps of an algorithm for 
determination of the degree of asymmetry of a defined DNA sequence using 
the calculation of DPT frequencies. The flow chart in figure 7 illustrates a 
further, even more specific example, of a preferred embodiment of the 
present invention, using the calculation of DPT frequencies, in the form of 
an algorithm implemented in PERL programming language. The variable 
names and functions indicated in bold in figure 7 are used by way of 
example and no details in these examples should be taken as limiting the 
application of this invention. The first step (101) is for a nucleotide sequence 
of a single strand of DNA (input sequence, $input_seq) of a desired length to 
be input. The sequence may be of any length from two bases to 3 X 10^ 
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bases, preferably from 5,000 to 50,000, and most preferably from 10,000 to 
20,000. The second step (102) is the input of length of the desired window 
size ($win_sym). Window size (W as described hereinabove) may be any 
number from one to a value equal to that of the largest whole integer smaller 
than one half the length of the DNA sequence, preferably from 20 to 300 and 
most preferably from 80 to 100. Beginning then at 103, the algorithm starts 
at the first potential axis of symmetry (axis position = 2* the window size) 
and calculates and outputs a set of DPT frequencies for that axis (104). The 
algorithm then calculates and outputs DPT residuals and one chi-square 
value per axis as described hereinabove. As for figure 4 as described 
hereinabove, the input sequence is converted to the two complementary 
sequence indexed arrays: @trgt_fwd and @trgt_revcomp in the steps 
indicated as 103 using string $win_seq in the process of these steps. In the 
next steps 105 and 105', the algorithm advances to the next potential axis; 
variables $basefeed and $basefeed_comp are used to advance the axis and 
sequence in the example of figure 7. In steps 106 and 106', DPT frequencies, 
residuals and chi-square values around the new axis are calculated and 
output and steps 105/105' and 106/106' are repeated iteratively up to and 
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including axis position = 2*length of the input sequence -(2*window size). 
In step 107 the ordered arrays of DPT frequencies, DPT residuals and 
chi-square values around each potential axis of symmetry are saved to a file. 
Further statistical, step 108, and graphical (step 109) analysis can be 
performed, allowing for identification of functional elements, step 110. 

Further envisioned as being within the scope of the present invention 
are alternate configurations wherein the nucleotide sequence of a specified 
window size upstream of the tested potential axis of dyad symmetry is 
compared against the nucleotide sequence of an equal size downstream firom 
this axis on the same, rather than the complementary, DNA strand. In the 
preferred embodiments illustrated in Figs. 3 and 6, a dyad pair on 
complementary strands is analyzed in order to determine s, and therefore p, 
or in order to determine DPT frequencies. Alternatively the sequence of only 
one single strand can be analyzed, based on the complementary nature of the 
strands. As a non-limiting example, rather than checking to see if the two 
bases that are situated at the same position relative to a potential axis of 
symmetry, but located on opposite strands of double stranded DNA, are 
identical, in this configuration the two bases that are situated at the same 
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position relative to a potential axis of symmetry, but located on the same 
DNA strand, (referred to as mirror pair), are examined to check whether they 
are complementary. For example, a ^ - g dyad pair and a ^ - ^ mirror pair 
are the same entity. 

Local biases in nucleotide composition can strongly contribute to dyad 
pair type frequency deviation because the frequencies of dyad pair types are 
proportional to the frequencies of occurrence of the bases from which they 
are comprised. For example, a GC rich region of DNA will have higher 
frequencies of the ° - g and ° - c dyad pair types. Thus, specifically 
envisioned as being within the scope of the present invention are alternate 
configurations wherein the frequencies of bases within a given region of 
DNA sequence of a defined window size on either side of a potential axis is 
determined. Such a method of determining nucleotide composition 
frequencies is inferior in accuracy to directly determining DPT frequencies 
because it neglects the effect of base order and therefore captures less of the 
available information than the direct DPT frequency analysis. Nucleotide 
composition frequency analysis can however be used for the same purpose 
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of locating functional elements and evaluating electrical conductivity, in a 
very similar way to the method described, albeit in a less definitive manner. 

Figure 10 is a high level block diagram of a system 30 for predicting 
the electrical conductivity properties and for identifying functional elements 
in a defined DNA sequence according to the present invention. System 30 
includes a processor 32, a random access memory 34 and a set of 
input/output devices, such as a keyboard, a floppy disk drive, a printer and a 
video monitor, represented by I/O block 36. Memory 34 includes an 
instruction storage area 38 and a data storage area 40. Within instruction 
storage area 38 is a software module 42 including a set of instructions which, 
when executed by processor 32, enable processor 32 to calculate dyad pair 
type frequencies, perform statistical analyses and graphical plotting by the 
method of the present invention. 

Using the appropriate input device 36 (typically a floppy disk drive), 
source code of software module 42, in a suitable high level language, for 
calculating dyad pair type frequencies, and performing statistical analyses 
according to the present invention is loaded into instruction storage area 38. 
The source code of software module 42 is provided on a suitable computer 
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readable storage medium 44, such as a floppy disk or a compact disk. This 
source code is coded in a suitable high-level language. 

Selecting a suitable language for the instructions of software module 
32 is easily done by one ordinarily skilled in the art. The language selected 
should be compatible with the hardware of system 30, including processor 
32, and with the operating system of system 30. If a compiled language is 
selected, a suitable compiler is loaded into instruction storage area 38. 
Following the instructions of the compiler, processor 32 turns the source 
code into machine-language instructions, which also are stored in instruction 
storage area 38 and which also constitute a portion of software module 42. 
Using the appropriate input device 36, the parameters of the DNA sequence 
analysis are entered, and are stored in data storage area 40. The results of the 
analysis are displayed at video monitor 36 or printed on printer 36. 

While reducing the present invention to practice, several results 
pertaining to the identification of functional elements in DNA sequence were 
obtained. These results clearly show that the method and system of the 
present invention identifies functional elements such as transcription related 
functional elements including promoters and enhancers and identifies 
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alterations related to regulatory domain mutations underlying human genetic 
diseases. 

Additional objects, advantages, and novel features of the present 
invention will become apparent to one ordinarily skilled in the art upon 
examination of the following examples, which are not intended to be 
limiting. 

Graphs of DPT frequency variation in several human genomic 
fragments, each containing a disease- associated gene, are presented in Fig. 8. 
Shown are output plots of dyad pair type frequency analyses of nine human 
genomic sequence fragments each containing a condition associated gene 
and one control fragment of computer generated random sequence. Each 
fragment is 40 kilobases (kb) long. A window of length 300 basepairs (bp) 
was used in all analyses shown. The start site of the primary transcript of 
each gene is marked by a yellow circle on the x-axis of the graph to which it 
corresponds, with an arrow indicating the direction of the gene. The x-axis 
represent bp positions on the + strand of the GenBank entry. The horizontal 
yellow line on each graph corresponds to zero on the left Y scale (number of 
residual dyad pairs); the black line corresponds to 30.6 on the right Y scale. 
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which is the value of the threshold x (15 degrees of freedom.; a=0.01). At 
residual values of zero, there is no difference between the observed and 

expected DPT frequencies; above, there is overexpression of that DPT, and 

2 

below, there is a supression of that DPT. All axes with x values above the 
30.6 threshold have, around that axis, DPT frequencies that deviate 
significantly from the expected. The sequences presented are: FMRl: 
Fragile-X mental retardation (GenBank accession #L_29074 , bp l-40k 
analyzed); WRN: Werner Syndrome (GenBank accession #181896, bp l-40k 
analyzed); POU4F3: Hearing impairment (GenBank accession #NT_006700, 
bp 120k- 160k analyzed); ATM: Ataxia Telangiectasia (GenBank accession 
#U82828, bp l-40k analyzed); RB: Retinoblastoma (GenBank accession 
#L11910, bp l-40k analyzed); NPCl: Niemann Pick CI syndrome 
(GenBank accession #NT_011044, bp 220k-260k analyzed); CFTR: Cystic 
Fibrosis (GenBank accession #AC_0001 1 1, bp l-40k analyzed); HEXA: Tay 
Sachs syndrome (GenBank accession #NT_0 10303, bp 190k-230k 
analyzed); and HD: Huntington disease (GenBank accession #88756, bp 
l-40k analyzed). 
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Some important features, common to all the graphs in Fig 8, are: (a) 
Significant deviation of DPT frequencies exists in the majority of axes in the 
sequences presented (all points with x values above the a=0.01 threshold 
line are significant), (b) In all of the sequences shown, the ^— a and ^—j 
DPTs are over represented and the c and ^— q DPTs are suppressed, and 
(c) There appears to be a defined pattern of DPT frequency variation, which 
coincides with the transcription start sites of genes. A dashed line circle is 
drawn around the DPT frequency variation element, coinciding with the 
transcription initiation site of each of the genes shown. This DPT variation 
element consists of a local, steep increase in the frequencies of the °— c 
and/or g DPTs, with a concomitant decrease in the frequencies of the ^— a 
and T DPTs. The length of the element is typically 1-2 kb, but is widely 
variable from gene to gene. 

Over 50 genes were analyzed so far; in the overwhelming majority of 
them a recognizable variant of the element was found to reside on the 
transcription start site. Although variations in the element are considerable 
from gene to gene, it appears that the main features of the element are 



49 

conserved. In addition to the transcription start site element, there are other, 
distinct, elements distributed in genomes of various organisms including 
viruses, bacteria and eukaryotes. DPT frequency deviation is indicative of 
anisotropy in an underlying physical property of the double helix, a property 
5 that is related to charge conductivity. Regions of DNA that possess DPT 
frequency deviation thus are used to identify both regions with altered 
electrical conductivity as well as functional elements within the DNA 
sequence. 

The present invention further permits identification and understanding 
10 of mutations underlying genetic conditions. Fig. 9 shows a DPT frequency 
analysis of two 40 kb sequences, with a window length of 300 bp. In the 
lower graph 9B, the scales of the Y axes are identical to those in the upper 
graph 9A, although the maximal values of chi square in graph 9B far exceed 
the maximal value on the axis and reach their maximum at >1800. The 
15 analysis in graph 9A is of the FMRl fragment containing bases 235k-275k 
from the GenBank entry accession #NT_0 11744. In graph 9B, the DPT 
analysis is of the "FMR1+(CGG)333" fragment which was obtained by 
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inserting a 1000 bp fragment (containing 333 tandem repeats of the (CCG) 
trinucleotide) into the sequence described for graph 9A, at position 255459 
of the GenBank entry #NT_01 1744. The fragment was inserted at the site of 
the (COG) repeat, expansion of which was shown to cause Fragile-X 
Syndrome. The "FMR1+(CGG)333" fragment thus simulates an expanded 
allele with approximately 350 (CGG) repeats. When affected and 
non-affected alleles of the gene known to be responsible for Fragile-X 
syndrome (FMRl) are analyzed and compared using the present invention, a 
clearly visible difference is found between the transcription initiation site 
elements of the affected and non-affected alleles. As many genetic conditions 
may be caused by mutations in regulatory control regions, the system and 
method of the present invention will not only permit identification of such 
mutations but will also therefore contribute to a better understanding of the 
control of gene expression in both normal and pathological situations. 

Although the invention has been described in conjunction with 
specific embodiments thereof, it is evident that many alternatives, 
modifications and variations will be apparent to those skilled in the art. 
Accordingly, it is intended to embrace all such alternatives, modifications 
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and variations that fall within the spirit and broad scope of the appended 
claims. 
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